Loss of coherence in dynamical networks: spatial chaos and chimera states 
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We discuss the breakdown of spatial coherence in networks of coupled oscillators with nonlocal in- 
teraction. By systematically analyzing the dependence of the spatio-temporal dynamics on the range 
and strength of coupling, we uncover a dynamical bifurcation scenario for the coherence-incoherence 
transition which starts with the appearance of narrow layers of incoherence occupying eventually 
the whole space. Our findings for coupled chaotic and periodic maps as well as for time-continuous 
Rossler systems reveal that intermediate, partially coherent states represent characteristic spatio- 
temporal patterns at the transition from coherence to incoherence. 
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Understanding the dynamics on networks is at the 
heart of modern nonlinear science and has a wide ap- 
plicability to various fields [J [5]. Thus, network science 
is a vibrant, interdisciplinary research area with strong 
connections to physics. For example, concepts of theoret- 
ical physics like the Turing instability, which is a known 
paradigm of non-equilibrium self-organization in space- 
continuous systems, have recently been transferred to 
complex networks [3J. While spatially extended systems 
show pattern formation mediated by diffusion, i.e., local 
interactions, a network takes also into account long-range 
and global interactions yielding more realistic spatial geo- 
metries. 

Network topologies like all-to-all coupling of, for in- 
stance, phase oscillators (Kuramoto model) or chaotic 
maps (Kaneko model) were intensively studied , and nu- 
merous characteristic regimes were found [IHl]. In par- 
ticular, for globally coupled chaotic maps they range - 
for decreasing coupling strength - from complete chaotic 
synchronization via clustering and chaotic itineracy to 
complete desynchronization. The opposite case, i.e., 
nearest-neighbor coupling, is known as lattice dynam- 
ical systems of time-continuous oscillators, or coupled 
map lattices if the oscillator dynamics is discrete in time. 
These kinds of networks arise naturally as discrete ap- 
proximation of systems with diffusion and have also been 
thoroughly studied. They can demonstrate rich dynamics 
such as solitons, kinks, etc. up to fully developed spatio- 
temporal chaos [5l410j. 

The case of networks with nonlocal coupling, however, 
has been much less studied in spite of numerous applica- 
tions in different fields. Characteristic examples pertain 
to neuroscience [HKT^, chemical oscillators [131 [H], elec- 
trochemical systems [IS], and Josephson junctions [IB] . 
A new impulse to study such networks was given, in 
particular, by the discovery of so-called chimera states 
[T71 [T5] . The main peculiarity of these spatio-temporal 



patterns is that they have a hybrid spatial structure, par- 
tially coherent and partially incoherent, which can de- 
velop in networks of identical oscillators without any sign 
of inhomogeneity. 

In this Letter we discuss the transition between coher- 
ent and incoherent dynamics in networks of nonlocally 
coupled oscillators. We start with coupled chaotic maps 

]=i-p 

where Zi are real dynamic variables {i = 1, TV, ^ 1 
and the index i is periodic mod N) , t denotes the discrete 
time, <T is the coupling strength, P specifies the number of 
neighbors in each direction coupled with the i-th element, 
and f{z) is a local one-dimensional map. We choose / as 
the logistic map /(z) = az(l — z) and fix the bifurcation 
parameter a at the value a — 3.8. This choice yields 
chaotic behavior of the map / with positive Lyapunov 
exponent A w 0.431. 

Results of direct numerical simulation of the model ([T]) 
in the two-parameter plane of the coupling radius r = 
P/N and coupling strength a are presented in Fig. [T] 
This figure reveals the appearance of regions of spatial 
coherence, shown in shading (color), at an intermediate 
radius of coupling. Alternatively, if the oscillators are un- 
coupled (r = 0) or coupling is only local (r = 1/N) the 
network displays high-dimensional space-time chaos [5]. 
In the opposite situation, when the coupling is all-to-all 
(r = 0.5), chaotic synchronization occurs: the oscilla- 
tors behave identically, but chaotically in time following 
the dynamics of / [3 [3]. The chaotic synchronization 
(hatched region fc = 0) persists for smaller r or tr up 
to the blowout bifurcation [TH] indicated by the curve 
BB, where the synchronized state loses transverse sta- 
bility, i.e., the dynamics becomes desynchronized. The 
spatially homogeneous state represents the simplest ex- 
ample of coherent dynamics. 



2 




0.2 



0.3 



A. /\ 



0.3 



i 100 i 100 



.A .A 


1 






n" 












.'■ 0.3 





i 100 



0.1 



0.2 0.3 0.4 0.5 

r 



FIG. 1: (Color online) Regions of coherence for system ([TJ in 
the {r,a) parameter plane with wave numbers k — 1,2, and 
3. Snapshots of typical coherent states Zi are shown in the 
insets. Color code inside the regions distinguishes different 
time-periods of the states. The CIB bifurcation curve sepa- 
rates regions with coherent and incoherent dynamics. In the 
hatched and shaded (color) regions below CIB two-cluster in- 
coherent states exist. Completely synchronized chaotic states 
exist in the light hatched region bounded by the blowout bi- 
furcation curve BE. Parameters: a = 3.8 and A'^ — 100. 



In general we call a network state zj, i = I, ...,iV, co- 
herent on the ring 5^ as iV — ^ c» if for any point x ^ 
the limiting value 



lim lim 

Af— i-oo t— i-oo 



sup 



z] I ^ 0, 



for (5-^0, (2) 



where Uf {x) = {j : < j < TV, \jlN - x\ < 8} denotes 
a network-neighborhood of the point x. If the limit ([2| 
does not vanish for (5 — 0, at least for one point a;, the 
network state is considered incoherent. 

Geometrically, coherence means that in the thermody- 
namic limit TV — > oo snapshots of the state z\ approach a 
smooth profile z(a;, i) of the spatially continuous version 
of Eq. (0 given by 



z'+\x)^!{z\x)) + 



2r 



[f{z'{y))~f{z\x))]dy. 

__ — r 

, (3) 

According to the definition given above, a transition from 
coherence to incoherence occurs when the respective so- 
lution profile z*(a:) becomes discontinuous. Note that for 
networks of phase oscillators the property of coherence 
and incoherence can also be established with use of the 
notion of a local order parameter ^20l |2T] . 

Regions of coherence and typical shapes of the respec- 
tive coherent states z* are shown in Fig. [l] as shaded 
(color) tongues and insets, respectively. A coherent state 
has a smooth profile characterized by the number of max- 
ima, i.e., the wave number k. Only regions for wave 
numbers k = 1,2, and 3 are shown. Further decrease 
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FIG. 2: (Color online) Coherence-incoherence bifurcation 
for coupled chaotic logistic maps for fixed coupling radius 
r — 0.32 (black triangles in Fig. [T]). For each value of the 
coupling parameter cr (decreasing from top to bottom, a = 
0.43,0.4,0.32,0.3,0.2, and 0.1, respectively) snapshots (left 
columns) and space-time plots (right columns) are shown. 
Other parameters as in Fig. [T] 



of r yields additional thin higher-order regions follow- 
ing a period-adding cascade k — 4,5,... Inside the re- 
gions, the states are coherent in space and periodic in 
time and undergo a period-doubling cascade of bifurca- 
tions in time as r or cr decrease. In the parameter space 
between the coherence regions the network dynamics re- 
main coherent but not periodic anymore. The states 
alternate chaotically between the adjacent fc-states and 
thus exhibit chaotic itineracy [51 [55] . The combination of 
period- adding in space and period- doubling in time repre- 
sents a remarkable feature of networks of coupled chaotic 
oscillators with nonlocal coupling. 

A typical scenario of the coherence-incoherence transi- 
tion is illustrated in Fig. [2]^ a)- (f), where we fix the cou- 
pling radius r — 0.32 and decrease the coupling strength 
a along the vertical line with triangles in Fig. [l] First, 
in Fig. [2|a) , the solution profile z* is clearly smooth for 
a = 0.43. Thus, the network dynamics is spatially co- 
herent. For smaller a, the profile z* sharpens up and. 
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at some value a = 0.40, loses smoothness in two points 
xi and X2 as shown in Fig. ^h). This is a bifurcation 
point for the coherence-incoherence transition: Beyond 
this parameter value, the wave-like profile z| splits up 
into upper and lower branches, and two narrow boundary 
layers of incoherence are born around the points xi and 
X2 (shaded yellow stripes ai and a2 in Fig. [2]jc)). The 
incoherence stripes become wider with further decrease 
of <7 (Fig. [2](d)) and, eventually, the dynamics becomes 
completely incoherent (Figs. [2je) and (f)). 

In our numerical simulations, no coherent states were 
found below the bifurcation parameter value a = 0.40. 
In contrast, numerous hybrid states arise, which are co- 
herent on some intervals of the ring and incoherent 
on the complementary intervals. Typical examples of 
these partially coherent states are shown in Figs. [2jc) 
and (d). In the figure, black diamonds mark a thresh- 
old within the incoherent regions: If the initial value for 
a chosen oscillator is located above/below this diamond, 
with all other oscillators unchanged, it will be attracted 
by the upper/lower branch. This implies that within the 
incoherent intervals ai and a2 any combinations of the 
upper and lower states- so-called mosaic [8 or skeleton 
pattern [10] - are admissible and can be realized by ap- 
propriate choice of the initial conditions. The coherence- 
incoherence transition is illustrated by the local order pa- 
rameter Ri shown in Fig. [sj^a) for the snapshots depicted 
in Fig. [2] It is defined as (cf. Ref. [H]) 



= lim 
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{z = l,...,N) (4) 



with the phase ipj introduced by the mapping sintpj = 
{2zj — maxj Zj— min^ Zj)/{'niaxj Zj—mhij Zj) and 5/N 
for — cx) , such that a spatial half-period oscillation is 
mapped onto the polar angular interval [— 7r/2, 7r/2]. Ri 
is close to unity for the coherent state and decreases in 
regions of spatial incoherence. For the two cases of com- 
plete incoherence (Fig.[2je),(f)) the local order parameter 
is much smaller than unity, and fluctuating strongly as a 
signature of spatial chaos. In case of very small coupling 
(Fig.[2];f)) the values of Zi are more spread out, and hence 
Ri varies less strongly between neighboring sites and is 
on average larger than in Fig.[2]je). 

As it is illustrated in the space-time plots of Fig. [2] 
the temporal dynamics before and after the coherence- 
incoherence bifurcation remains periodic up to very small 
coupling, when it is chaotic (Fig. [2jf)). The system's 
complexity results from the fact that the bifurcation gives 
rise to a huge multistability of partially coherent states 
as — ?> oo. Indeed, the number cjv of different par- 
tially coherent states born in the bifurcation is cat = 2''^, 
where d is the fraction of oscillators in the incoherent part 
[d = ai +Q!2 in the case of two incoherence intervals as in 
Figs.[2]jc) and (d)). It follows that the number of different 
states grows exponentially fast with N , and the spatial 
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FIG. 3: (Color online) (a): Local order parameter for r = 0.32 
and different g as in Fig. [2] (approximated with 5 = IQ, N = 
400). (b): Regions of stability for the two-cluster solutions 
with different asymmetries (ni : 712) as a function of o for 
r = 0.35. Other parameters as in Fig. [l] 



entropy h, which is defined as h = lim7v-i.oo(l/-^) In c^r, 
is positive and equals h — din 2. Positive spatial entropy 
means that the system displays spatial chaos [5HlO). i.e., 
sensitive dependence on space coordinates. Therefore, 
the coherence- incoherence bifurcation results instantly in 
the appearance of spatial chaos that develops first at nar- 
row incoherence intervals and, with decreasing a, spreads 
onto the whole ring. Thus a chimera-like state of coex- 
isting coherent and incoherent regions arises as a tran- 
sitional state in the coherence-incoherence bifurcation 
scenario. However, in contrast to previously reported 
chimera states in time-continuous systems [171 118] , the 
temporal behavior is periodic rather than chaotic, and 
the complexity arises due to the huge variety of multi- 
stable incoherent states corresponding to permutations 
of the sequence of upper and lower local states. With 
further decrease of cr, the chimera states disappear giv- 
ing rise to completely incoherent behavior. 

To identify the parameter range for partially coher- 
ent states, we define a mosaic of z* as a symbolic se- 
quence of "— " or if the value Zi belongs to the 
lower or upper branch of the solution profile, respectively. 
Hence, states as in Fig. [2|^b) are given by the mosaic of 

the form (... h +... + H ...). Therefore, they 

may be considered as two-cluster states with the ratio of 
(ni : 77.2), where ni and n2 are the numbers of "— " and 
(ni +n2 = N), respectively. The ratio (rii : ^2) in- 
dicates the level of asymmetry of the solution z* , which 
is important for its stability. Indeed, as it is illustrated 
in Fig.jsf^b) for r — 0.35, the symmetric solution (50 : 50) 
has the widest stability interval a G (0.189, 0.41). As the 
asymmetry grows, the stability interval shrinks, and the 
solution with the mosaic ratio (35 : 65) has the shortest 
stability interval cr £ (0.314,0.317). Two-cluster solu- 
tions with larger asymmetry cannot be stabilized any- 
more. The states with more complex mosaics, examples 
of which are presented in Fig. [2]^c)-(e), are characterized 
by more involved mechanisms of stability. 

To test if the coherence-incoherence bifurcation is a 
universal scenario, we have also investigated nonlocally 
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FIG. 4: (Color online) Snapshots with decreasing coupling 
strength a: Coherence-incoherence transition for a ring of (a) 
100 coupled periodic logistic maps (a = 3.2) for a coupling 
radius r = 0.1 and (b) 100 nonlocally coupled Rossler systems 
(a = 0.42, 6 = 2, c = 4) with r = 0.3. 



coupled networks with different local dynamics. Figure |4] 
shows the coherence-incoherence transition for nonlocally 
coupled logistic maps ([ij in the periodic regime (a = 3.2, 
period 2, Fig. |4][a)) and for nonlocally coupled chaotic 
Rossler systems 



Zi + 



+P 



2-P j=i~p 



Xi) . 



(5) 



yi = Xi+ ayi, 
Zi ^ h Zi{xi ~ c) (j = l,...,iV) 

(Fig. |4]jb)). As it can be seen from the snapshots, both 
models display a transition from spatial coherence to in- 
coherence, as the coupling strength a decreases, accord- 
ing to the scenario described above. Since network ([5| 
is time-continuous, the oscillators within the incoherence 
intervals are not only located at an upper or lower branch 
of the solution profile, but vary continuously. This gives 
rise to chaotic temporal dynamics in the incoherent inter- 
vals, which resembles known chimera states [TTIIIH]- We 
conclude that chaotic chimera states typically arise in the 
nonlocally coupled Rossler systems, similar to nonlocally 
coupled Kuramoto-Sakaguchi phase oscillators [18J. 

In conclusion, we have identified a novel mechanism 
for the coherence- incoherence transition in networks with 
nonlocal coupling of variable range. It consists in the 
appearance of multistable chimera-like states. We have 
found similar bifurcation scenarios for coupled maps with 
both chaotic and periodic local dynamics as well as 
for time-continuous systems. This indicates a common, 
probably universal phenomenon in networks of very dif- 
ferent nature, due to nonlocal coupling. 
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